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ABSTRACT 

Context. The multi-purpose INTEGRAL mission is continuously contributing to Gamma Ray Burst (GRB) science, 
thanks to the performances of its two main instruments, IBIS and SPI, operating in the hard X-ray / soft 7-ray domain. 
Aims. We investigate the possibilities offered to the study of GRBs by PICsIT, the high-energy detector of the IBIS 
instrument. 

Methods. We searched for transient episodes in the PICsIT light curves archive from May 2006 to August 2009, using 
stringent criteria optimized for the detection of long events. In the time interval under examination PICsIT provides 
an energy coverage from 208 to 2600 keV, resolved in eight energy channels, combined with a fine time resolution of 
16 ms. 

Results. PICsIT successfully observes GRBs in the 260 — 2600 keV energy range with an incoming direction spread over 
half the sky for the brightest events. We compiled a list of 39 bursts, most of which are confirmed GRBs or simultaneous 
to triggers from other satellites/instruments. We produced light curves with a time sampling down to 0.25 s in three 
energy intervals for all events. Because an adequate response matrix is not yet available for the PICsIT burst sample, 
we obtained a calibration coefficient in three selected energy bands by comparing instrumental counts with physical 
fluences inferred from observations with different satellites. The good time resolution provided by the PICsIT data 
allows a spectral variability study of our sample through the hardness ratio. 

Key words. Gamma-ray burst: general - Instrumentation: detectors 
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10 1- Introduction (|Rau et al.ll2005f) . detected with the anticoincidence shield 

(ACS) of the SPI instrument. These light curves, with 50 

', The INTEGRAL satellite (jWinkler et al.l [20031 ) . launched ms resolution, are used for GRB trianRulatio ns with other 

Q . in October 2002 and still in operation, is an interna- satellites of the inter planetary network (IPN, iHurlev et al.l 

T— I ' tional observatory exploiting a set of coded mask instru- 2004). Although SPI- ACS pr ovides no spectral information, 

T— I I ments to explore the sky in the hard X-ray / soft 7-ray IVigano fc Mereehetti] (|2009ll derived a calibration factor to 

^ ' domain with fine imaging and high spectral resolution, convert the instrumental count rate to physical flux on the 

] Although INTEGRAL was conceived as a multi-purpose basis of the observation of a large sample of GRBs. 

• ^ high-energy mission not specifically optimized for the ob- In this paper we expl ore the possibilities for GRB stud- 

^ ' servation of GRBs, it has been making a relevant contribu- ies offered by PICsIT ()Di Cocco et al.l [20031) . the high- 

H I tion to the study of these fascinating and in many respects energy detector of the IBIS instrument, by presenting the 

1 still mysterious events thanks to the high sensitivity and results of a search for GRBs detected from May 2006 to 

good imaging performanc es of its two main instrument s August 2009. 
IBIS (jUbertini et al.ll2003l) and SPI ([Vedrenne et al.ll2003l) . 
Gamma ray bursts are detected in the IBIS field of view at 

a rate of ^10 per year. Their positions, with an accuracy of 2. Data analysis 

few arcminutes, are automatically distributed in near real 2 j PICsIT data 
time thanks to the IN TEGRAL burst alert system (IBAS, 

iMereghetti et al.l 120031 ). enabling rapid follow-up observa- INTEGRAL moves along a highly eccentric orbit with a 

tions. A systematic spectral analysis of all GRBs detected revolution period around the Earth of 3 days, perigee of 

with I BIS un til 2008 has been p resented bv lVianello et al.l ~ 10,000 km and apogee of ~ 150,000 km. Observations 

([20091 ) ■ while [pblev et al] (|2008f ) carried out a comprehen- are carried out when the satellite is above the radiation 

sive study of a sample of 46 bursts, focusing mainly on belts, at an altitude higher than ~ 40, 000 km. To optimize 

their timing properties using IBIS and SPI. INTEGRAL is the image reconstruction a dithering strategy is adopted 

also providing light curves for about 15 triggers per month during most scientific acquisitions: each observation is split 
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Fig. 1. Count distribution and power spectrum in four energy bands for ScW 04820018. Count distributions are obtained 
with a time resolution of ~ 16 ms {first row), 1 s {second row) and 10 s {third row). Poissonian {upper panels) and Gaussian 
{second and third rows) models are superimposed as (red) dashed lines. The number of bins is 222256, 3462, 345, for 
16 ms, 1 s and 10 s, respectively. The power spectra {bottom panels) show a white noise background in all energy bands 
with the exception of the 208-260 keV channel. 



into a set of pointings (science windows - ScWs) covering a 
grid of directions spaced 2° apart according to a rectangular 
or an hexagonal pattern. The typical duration of each ScW 
is ~ 2, 000 s. 

The PICsIT detector consists of a 64 x 64 pixels array 
of CsI(Tl) scintillators, operating in the nominal energy 
band 175 keV-10 MeV. A combination of active and pas- 
sive shieldings partially surrounds the detector to reduce 



the background and to limit its field of view to the direc- 
tions corresponding to the coded mask (29° x 29°). With 
increasing energy the passive shielding becomes progres- 
sively transparent to the incoming radiation and PICsIT 
is expected to detect GRBs even at far off-axis directions, 
outside its nominal field of view. 

The first GRB detected by PICsIT (GRB 021125- 
iMalaguti et al.l [20031 ) occurred during the in-flight calibra- 
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Channel Energy Band 



(keV) 

1 208-260 

2 260-312 

3 312-364 

4 364-468 

5 468-572 

6 572-780 

7 780-1196 

8 1196-2600 



Table 1. On-board configuration of the energy channels 
of spectral timing data since May 2006. The time binning 
is 16 ms for all energy channels. 



tion phase, when data acquisition (in photon-by-photon 
mode) maintained the full spatial, temporal and energy 
information. When PICsIT is operated in standard scien- 
tific mode, its data are pre-processed on-board and orga- 
nized in different kinds of histograms before they are sent 
to the ground. This is required to comply with the tight 
telemetry rate available, which does not permit to trans- 
mit the whole information (time, energy and position) of 
each detected count. The standard histograms consist of 
spectral imaging (SI) and spectral timing (ST) data, op- 
timized for spatial and temporal resolution, respectively. 
Spectral imaging histograms preserve the full spatial infor- 
mation and can be used to derive sky images in different 
energy bands. However, they are integrated over the whole 
duration of each ScW (typically a few thousand seconds), 
and are in general useless for GRBs, which are bright on 
shorter time scales and therefore diluted and not detectable 
in the SI data (the only exception to date is the very bright 
and long GRB 04121 9A, which has been d etected in both 
SI and ST data sets - iDi Cocco eraI1l2007D . 

In this work we used the ST data, which provide the 
total rate of events on the detector plane with a configurable 
time binning down to 1 ms, and an energy resolution in 
up to eight energy channels. The ST data do not permit 
to produce images because the spatial information is not 
preserved. The on-board configuration of the ST data has 
been changed a few times during the mission lifetime to 
refine the energy sampling and enhance the signal-to-noise 
level of the detected sources. We restricted the analysis to 
the time period from May 2006 to August 2009 (revolutions 
440-834), during which the ST time binning was set to 
16 ms and the energy channels, covering the range 208- 
2600 keV, remained fixed at the nominal values given in 
Table [3 

2.2. GRB search 

In order to search for possible GRBs detected in the PICsIT 
ST data, we examined the light curves looking for statisti- 
cally significant excesses above the background counts. The 
detector background level is not constant: intra-revolution 
variations are visible when the satellite is close to the ra- 
diation belt passage, the relative position of strong sources 
can cause different background level as a function of the 
pointing direction, and solar activity can also produce back- 
ground variations. In addition, the background count rate 
shows a long-term evolution, which can be ascribed to the 



detector activation and to the variable cosmic ray environ- 
ment, modulated with the solar cycle. 

We performed our GRB search analyzing each ScW in- 
dividually, because these are the longest homogeneous pe- 
riods with a fixed satellite attitude, and are sufficiently 
short, compared to the timescales of orbital background 
variations, to approximate the background with a constant 
value. Therefore, we evaluated the expected background 
level as follows for any given ScW, energy channel and time 
rebinning. We first removed all time intervals with large 
spikes or drops verifying: \Ci— < C > \ > 6 ■ ac , where Ci 
is the number of counts in bin i and < C > and ac are the 
mean and standard deviation of counts, respectively, com- 
puted over the whole ScW. We then defined < B > and as 
as the mean and standard deviation over the cleaned ScW. 
Once the energy band and the time binning were selected, 
we defined a burst candidate as an excess in the light curve 
where the number of counts Ci is higher than a thresh- 
old given by < S > +Nth ■ ctb- The value of Nth and the 
sample size (in units of bins) establish the probability P of 
detecting statistical fluctuations above the threshold. 

In order to select the most convenient energy range and 
time binning for the GRB search, we inspected the distri- 
bution of counts within ScWs where no GRBs occur for 
different time binnings and in all available energy channels. 
Some examples of these count distributions are shown in 
Fig. [T] We found that the count distribution strongly de- 
viates from the Poissonian statistics in the lowest energy 
channel. This is confirmed by the Fourier analysis of the 
light curves. As shown in Fig. [1] (bottom panels), the power 
spectra above 0.2 Hz are consistent with a constant model 
(c = 2.23 ± 0.08, 2.31 ± 0.09, 2.17 ± 0.07, in 260 - 312, 
312 - 364 and 572 - 780 keV, respectively), with the ex- 
ception of the first energy band, which follows a power law 
trend above 0.1 Hz. Above 260 keV the count distribution 
is still not consistent with the expected Poissonian statis- 
tics, because power spectra describe a white noise process 
with observed variance exceeding the Poissonian of at least 
10%. For this reason the threshold level for the GRB search 
is expressed in terms of standard deviation. The different 
behavior of the first energy band can be ascribed to a sig- 
nificant contribution of charged particles interacting with 
the detector Csl crystals and inducing long-lasting phos- 
phorescence emissio n, which is not elimi nated by the anti- 
coincidence system (jSegreto et al.l [20031) . On the basis of 
this analysis we decided to perform the GRB search in the 
260 — 312 keV energy band, which seems to be the least 
affected by the systematic background component and, for 
the typical spectra of GRBs, maximizes the signal-to-noise 
ratio. 

A preliminary analysis was conducted on the data sam- 
ple to explore several rebinning and several threshold lev- 
els. We selected a threshold level Nth = 8 applied to light 
curves rebinned on 10 s time bins. This value guarantees 
a probability P < 10^^ of statistical fluctuation above the 
threshold, over the whole data set, of ~ 8 • 10^ time bins. 
The search with 10 s bins was carried out twice, with a 
half-bin offset. 

We searched for the GRB candidates using the pre- 
processed data of all ScWs, without an a priori selection to 
eliminate the ScWs with potential problems caused by, e.g., 
high solar activity, bright variable sources, high particle 
background close to the radiation belt entrance/exit, par- 
ticular instrumental settings or deviations from the nomi- 
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Table 2. Obscrvablcs of the PICsIT burst sample 



Burst 


ScW 


UT 


T90 


Fluence 


Peak Flux 


Energy range 


Instrument 


Ref 


Notes 






(hh:mm:ss) 


(s) 


(ct) 


(ct/s) 


(keV) 








060805B 


4650034 


14:27:15 


4.7± 0.4 


6315±186 


3600±208 


260-2600 


lA, K, R 


(1) 




060819 


4700017 


18:28:13 


14 ± 2 


12228±305 


2394±195 


260-2600 


lA, K 


(a) 


060901 


4740058 


18:43:56 


8 ± 1 


2199±182 


1671±186 


260-780 


II, K 


(2) 


(b) 


060905 


4750078 


14:48:51 


239 ± 47 


25055±849 


1918±188 


260-572 




060928 


4830025 


01:20:06 


31 ± 4 


38936±462 


4966±220 


260-2600 


B, lA, K, R, S 


(3) 


(c) 


061031 


4940040 


12:19:47 


34 ± 3 


4707±329 


1367±186 


260-572 


K 


(d) 


061122 


5010073 


07:56:52 


>1.50 


>2196 


> 3212 


260-780 


II, K 


(4) 


(e) 


061222A 


5110078 


03:30:15 


10.5± 1.5 


3509±236 


1871±191 


260-1196 


B, K 


(5) 


(f) 


070207 


5270055 


21:00:48 


lO.Oi 0.5 


6178±255 


3700±210 


260-2600 


B, lA, K, S 


(6) 




070227B 


5340043 


21:39:23 


32± 13 


3341±316 


672±179 


260-572 


lA, R 


(g) 


070326 


5430018 


00:45:26 


26.2± 0.7 


3080±302 


1466±189 


260-468 


lA 




(f) 


070329 


5440037 


14:59:48 


40± 3 


8720±495 


1048±185 


260-2600 


lA, K 




(g) 


070403 


5460003 


12:40:28 


47± 3 


5788±347 


701±181 


260-468 






070418 


5510013 


17:16:19 


31± 2 


5471±369 


1461±192 


260-780 


K 




(f) 


070429 


5540054 


02:40:28 


20.50±0.4 


3081±326 


1634±192 


260-2600 






(f) 


070829 


5950047 


20:08:36 


73± 1 


9766±618 


1762±194 


260-1196 


lA, K 




070917 


6010058 


04:41:24 


7± 2 


2946±202 


1119±95 


260-780 


lA, K 




(f) 


071003 


6070014 


07:40:54 


12.0± 0.3 


16175±302 


4013±218 


260-2600 


B, K, II 


(7) 


071006 


6080009 


06:42:08 


47± 8 


3332±401 


940±188 


260-572 


B, K, lA 


(8) 





Column description: Burst: name of the event (YYMMDD); ScW: pointing ID in which the event occurred; UT: trigger time; 
T90: duration of the event; Fluence: fluence in units of counts, referring to the Tgo interval and to the energy interval given in 
column 7; Peak Flux: in units of count/s in the same time and energy intervals as for fluence; Instrument: telescopes/instruments 
that observed the event (A = AGILE, B= Swift/BAT, F = Fermi/GBM, II = INTEGRAL/lSGRl, lA = INTEGRAL/SPl- ACS, 
K= Wind/ Konus, R = RH ESSI, S = Suzaku/WAM); R ef: GRB s have been documented by: (1) = iHurlev et all ll2006bl ). (2) 
= [M ereghet ti et al.l (l2006bl ). ( 3) = iHurlev et all mOQeli ) (4) = [Mereehetti et al.l (|2006al '). (5) = iGruoe et alJ (|2006l '). f6) = 
iGoienetsku et all (|20073 r77i = ISchadv et al.1 (|2007l ). (8) = ICummingsl (|2007l ). Notes: (a) = T90 refers to the 260-312 keV range, 
that gives a better signal to noise ratio, (b) = The bright GRB 060901, detected by INTEGRAL over ~20 s, saturated ISGRI 
telemetry but did not affect PICsIT data. In PICsIT light curve the burst shows a multipeak structure extending over 8 s, in 
accordance with Konus main event (followed, in Konus data, by a faint tail up to 15 s). (c) = GRB 060928 is structured in two 
main episodes at 01:17:05 and 01:17:20, with duration ~ 14 s and ~ 25 s, respectively. PICsIT observed both events (Fig. I13p but 
our analysis has been restricted to the main second event (Fig. I14|) . (d) — A linear fit of the background was required for T90 
co mputation, (e) = PICs IT light curves are severely affected by data gaps at all energies (Fig. I16p . A Tgo of 11 s was obtained 
bv lMcGlv nn et al.l (|2009l ') with SPI. (f) — For T90 computation the background was fitted to a polynomial function of the third 
degree, (g) = For T90 computation the background was fitted to a polynomial function of the second degree. 



nal operating parameters. The resulting ScWs with GRB 
candidates were then manually screened for the above con- 
ditions and re-analyzed with the s pecific standard analysis 
software osa v. 7.0, described in iGoldwurm et all ()2003[ ) 
to derive the relevant physical quantities of the events. The 
time of occurrence of each burst was refined using light 
curves with 1 s resolution and then cross-checked with solar 
fiares or high activity revealed by the INTEGRAL r adia- 
tion environment monitor (IREM, IHaidas et al. 20031). the 
GEOS neE and RHESSI monitoring (jLin et al.l 12001 to 
exclude particle events. 

We also performed a search with binning of 1 s and 
Nth — 12 on a subset of our data (revs. 440 — 680). This re- 
sulted in 21 events, 10 of which were also detected with the 
10 s bin search. However, a deeper analysis of the remain- 
ing 11 events suggests that they are not GRBs. Indeed, the 
light curves of all these events with the on-board resolu- 
tion (^ 16 ms) show a similar pattern, consisting of three 
peaks separated by 0.1 s, the central one dominating the 
high-energy bands and progressively dimming as energy de- 
creases (Figl2]) . Although the origin of these short events is 
uncertain, they are probably an instrumental effect caused 
by charged particles interactions. Since the search for short 
episodes did not add genuine GRBs to the sample of bursts, 
we restricted the search to the 10 s binning. 



200 ■ 

150 ■ 
100 ■ 
50 ■ 

88 
60 
I 40 
20 

108 
80 
60 
40 
20 
88 
60 

= 40 
J 

20 




05350023 



208-260 keV 




260-.31 2 kcV 



^^^^ 



_ 31 2-364 keV 



364-463 keV 




-1.0 -0.5 0.0 0.5 1.0 
Time since trigger (s) 



1 

^ 1 



. 468-572 keV 1 




572-780 keV 


L 


. 780-1 196 keV 








■ 1 196-2600 ke\ 








.0 -0.5 0.0 


0.5 1 



Time since trigger (s) 



Fig. 2. Typical light curve of short fake events, probably 
caused by instrumental effects. Light curves in eight en- 
ergy bands and 16 ms time resolution refer to the event in 
ScW 5350023. 



^ http://www.sec.noaa.gov 
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Table 2. Observables of the PICsIT burst sample - Continued. 



Burst 


ScW 


UT 


T90 


Fluence 


Peak Flux 


Energy range 


Instrument 


Ref 


Notes 






(hh:mm:ss) 


(s) 


(ct) 


(ct/s) 


(keV) 








071108 


6190025 


21:40:29 


20.0± 0.5 


9976±372 


5072±229 


260-2600 


lA, K, A 




(h) 


080122 


6440034 


18:32:43 


105± 4 


10796±743 


2112±199 


260-1196 


lA, K, S, A, B 


(9) 


080204 


6480067 


13:56:35 


14± 3 


9901±312 


4254±220 


260-2600 


K, lA, S, A, B 


(10) 




080303 


6580014 


21:34:45 


29. 2± 1.5 


31507±461 


8283±254 


260-2600 


lA 




080319B 


6630025 


06:12:50 


42± 1 


29392±536 


1905±198 


260-2600 


B, K 


(11) 




080328 


6660038 


08:03:12 


33± 2 


6593±419 


2243±200 


260-1196 


B, K, S, lA 


(12) 




080408 


6700004 


10:23:25 


40± 1 


6355±424 


1464±194 


260-780 






(i, f) 


080514B 


6820005 


09:55:58 


6± 1 


5611±195 


3495±214 


260-2600 


A, K, S,IA 


(13) 


080607 


6900002 


06:07:27 


10± 1 


4748±195 


2225±201 


260-2600 


B, A, K, lA 


(14) 




080613B 


6920007 


11:12:40 


31.7± 0.5 


5181±417 


1712±197 


260-1196 


B, K, lA 


(15) 




080615 


6920047 


04:07:33 


9.5± 0.3 


6769±239 


2768±206 


260-1196 


K, R 






080721 


7040083 


10:25:14 


16.0± 0.3 


14627±340 


4888±227 


260-2600 


B, K, lA, R 


(16) 




080723B 


7050036 


13:22:34 


45± 2 


>9936 


>3030 


260-780 


II, A, K 


(17) 


(j) 


080817A 


7130052 


03:52:16 


42± 2 


10461±528 


1509±196 


260-1196 


F, K, lA 


(18) 




080918 


7240065 


09:44:37 


2.5± 0.3 


3986±147 


3916±220 


260-2600 


K, lA 






090528B 


8080058 


12:22:57 


81 ± 7 


6665±622 


1141±199 


260-780 


F, S, K, lA 


(19) 


(f) 


090618A 


8150054 


08:29:26 


15± 7 


3032±301 


865±198 


260-1196 


B,A,F,K,S,IA 


(20) 


(k, f) 


090623 


8170026 


02:34:30 


40± 2 


6517±489 


1540±204 


260-1196 


F, S, K, lA 


(21) 




090626A 


8180029 


04:32:11 


42± 1 


13167±545 


3577±221 


260-2600 


F, S, K, lA, B 


(22) 




090809B 


8330036 


23:28:17 


9± 1 


3399±257 


1501±203 


260-1196 


F, K lA 


(23) 


(f) 



Ref ■ G RBs have been do cumented _byU^ = Hurlcv ct all (|2008D. (10) = LG olcnctsk ij_et all (l2008j'). fll) =LR acusin et al.' (|200i 'l 
(12) = iPerri et all ()2008), ( 13 ) = iRapisarda et a l. (200j) 7(14) = iMangano ct al. ( 20M). (151 = iMarkwa r dt et al.. (,2008'). (16) = 
[MarshaU ct al. (2008), (17) = lGotz et al.l ll2008fl. (18) = iBissaldi ct al. (2008). jl9) = fvon KienlinI (|2009aD . (20) = ISchadv et all 
(2009 ). (21) = ,Rau (,2009 ), (22) = lvon KienlinI l|2009bf ). (23) = van der Horst (,2()09l ). Notes: (h) = We analyzed the main pulse of 
the event, which is composed of three peaks over ~70 s. (i) = The burst was observed by PICsIT in two main episodes separated 
by 50 s and extending over 120 s. The first smaller peak seems to be affected by count drops and therefore the analysis was 
conducted on the second event only, (j) = The bright GRB 080723B fell in the ISGRI FoV, with an off-axis angle ~ 8°, and 
saturated the available telemetry. PICsIT data show gaps at all energies that are limited to few time bins, however, and did not 
affect the T90 computation. We report the lower limits of the peak flux and fluence. (k) = in PICsIT data the burst occurs ~60 s 
after the Swift/BAT trigger and coincides with the brightest portion of this long event (~ 130 s). 



3. Results 

The PICsIT burst sample selected with the 10 s binning 
as described above in the period May 2006 - August 2009 
consists of the 39 events listed in Table [2j Most of these 
events (23) have also been detected by other satellites, 
which provided their localizations, confirming their GRB 
nature. We are confident that the other events, for which 
no localizations have been obtained, are also most likely 
GRBs. Indeed, 11 of them have also been detected by other 
instruments such as the INTEGRAL/SPl-ACS (sensitive 
to photons above ^^80 kcV), Wind/Konus (energy range 
10 keV-10 MeV), and RHESSI (50 keV-17 MeV). Only four 
events were detected exclusively by PICsIT. Sixteen bursts 
are hard events detected up to the highest energy channel 
(1196-2600 keV), 26 reach the 780-1196 keV band, 5 are 
not visible above 572 keV, and 2 are observed only up to 
468 keV. 

We evaluated the background level in each energy chan- 
nel by fitting a polynomial to counts measured in two time 
intervals of a few hundred seconds manually selected be- 
fore and after the GRB. In all cases a constant was ad- 
equate, while a polynomial function up to the third de- 
gree was required for a better determination of the T90 
interval of 12 events (see notes in T able [21) ■ The Top du- 
ration was computed as described in iKoshut et al.l ( 19961 ) 
from the background-subtracted light curves in the energy 
range 260-2600 keV, with the exception of 060819, whose 
duration was evaluated in the 260-312 keV light curve. The 
background-subtracted light curves in three energy bands 
are reported in Figs. 



The values of fluence and peak fiux in Table [5] are re- 
ported in instrumental units. Indeed, in the next section we 
argue that a large part of bursts is expected to occur far 
off- axis, where a response matrix is not available yet. The 
fiuences refer to the Tgo intervals and to the energy range 
reported in column 7: the energy coverage for each event is 
given by the energy channels where the fluence F verifies 
the condition F > 2 ■ ap, where ap is the error on the flu- 
ence. The first energy channel (208-260 keV) was excluded 
because it is dominated by the systematic background com- 
ponent described in the previous section. The peak fluxes 
in the same energy range are for an integration time of 1 s. 



3.1. Spectral characterization and PICsIT effective area 

In Fig. [3] we show the cumulative distributions of the in- 
strumental fiuences and peak fiuxes of the PICsIT burst 
sample in three energy intervals: 260-364 keV, 364-780 keV 
and 780-2600 keV. The irregular behavior is likely caused 
by a non linear relation between instrumental observables 
and physical fiuences and peak fiuxes. 

Although not provided by the ST data, the arrival direc- 
tion of the photons can be expressed in PICsIT coordinates 
for GRBs with known sky positions from observations with 
other satellites. This can be made using the instrument atti- 
tude parameters stored in the housekeeping data on-board. 
In Fig. m we show the spatial distribution in PICsIT coor- 
dinates of 22 events in our sample. The polar angle Q is 
the off-axis angle, while is the azimuth, defined so that 
the direction toward the Sun is at = 90°. In the figure we 
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Fig. 3. Integral distribution of fluence (upper panel) and 
peak flux (lower panel) in instrumental units of the PICsIT 
burst sample. 



Fig. 5. Fluence 260-2600 keV as a function of instrumental 
coordinates Q (upper panel) and (j) (lower panel) for the 
21 GRBs with known spectral parameters . The telescope 
aperture is shown as a dotted line in the upper figure. 
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Fig. 4. Instrumental coordinates of 22 GRBs in the PICsIT 
sample. The (red) square gives the approximate PICsIT 
FoV and dashed line borders the area where PICsIT is par- 
tially shielded by the spectrometer SPI. 



indicated the IBIS FoV and the approximate region occu- 
pied by the spectrometer SPI (0 £ [250°, 290°], 6 > 28°). 
Despite the low number of events, we note that the region 
that includes the directions shielded by SPI in the {&, (j)) 
plane is scarcely populated. 

For 21 GRBs in our sample the fluences in physical units 
and the average spectral parameters are available from ob- 
servations with other satellites. The spectrum of six events 
in this subsample is described by a cut-off powerlaw, while 
15 events follow a Band function with mean parameters 
a = -0.9 ± 0.2, P = -2.5 ± 0.4, Ep = 318 ± 174 keV. On 
the basis of the spectral parameters reported for each sin- 
gle GRB, we converted the fluence to the PICsIT energy 
range (260-2600 keV), assuming no spectral evolution and 
neglecting differences in the reported burst duration. When 
the same event was measured by more than one instru- 
ment, we chose the one whose energy coverage overlapped 
the most with that of PICsIT. The derived fluences are re- 
ported in Table [3] and plotted in Fig. [S] as a function of 
PICsIT angular coordinates. We note that the four bursts 
within the PICsIT FoV (9 < 15°) span about one order 
of magnitude in fluence, while at wider off-axis angles only 
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Fig. 6. Ratio of fluences in physical and instrumental units 
in three energy bands for a subset of 20 GRBs. In the lower 
panel PICsIT counts were corrected for the geometrical fac- 
tor cos(0). 



Table 3. Position in instrumental coordinates and fluence 
in 260-2600 kcV for the subsample of documented GRBs. 





Burst 


e 





F (260-2600 keV) 
(10^^ erg cm"^) 


Ref. 






n 


(°) 






060805B 


72.4 


359.3 


5.9±0.2 


(1) 




060901 


14.1 


297.4 


76+°-?^ 

^' ' ^ — 0.10 


(2) 




060928 


62.4 


43.2 


iq 17+103 


('3') 




061122 


7 A 


990 Q 


r. C9+0.01 
"•^^-0.03 






061222 


22.3 


226.5 


1 K+0.2 


(5) 




071003 


13.5 


273.1 


3-8Io.5 


(6) 




071006 


55.6 


233.3 




(7) 




080122 


48.0±0.3 


240±1 


4±2 


(8) 




080204 


18±1 


101±2 




(9) 




080319B 


79.7 


215.7 


39.5lJ:° 


(10) 


080328 


22.6 


0.9 


l.OiO.l 


(11) 


080514B 


37.7 


355.4 


1.5±0.1 


(12) 


080607 


73.7 


201.4 


4.8±0.2 


(13) 




080613B 


o9.4 


io4.z 


i.oitU.o 


(14) 




080721 


58.2 


353.9 


5.0±0.4 


(15) 




080723B 


7.9 


104.7 


4.3±0.3 


(16) 




080817A 


57.7±0.6 


94±1 






090528B 


32±1 


215±1 


2.29±0.07 


(17) 




090618A 


85.4 


14.6 


9.8±0.2 


(18) 




090623 


47±1 


244±1 


^•"-0.9 


(19) 




090626A 


46±1 


182±1 


3.1±0.4 


(20) 




090809B 


53.7 


41.1 


1.63±0.02 


(21) 



The errors in polar coordinates were neglected when 
< 0.1°. Fluences were extrapolated in the 260-2600 keV 
energy rang e from spectra l parameters rep o rted in : (1) — 
BcUmctal. (200^, (2) = iGolenetskii "eraT (2006a'), (3) = 
Golcnctskii ct al. (2006y|), (4) =[Gplcnctskii ct al. (2006c), (5) = 
IGolen etskii et al. (2006d), (6) = lGolen ctskii ct al. (2007a), (7) 
IGolenetskii et al. (2 007b'). (8) = IGolenetskii et al. (20Q8ai), 
iGolenetskirit a l. (2008^), (10) = IGolenetskii etal.l 



iGolei 

= IGolenetskii et al. (2 007bl. (8) = 
(9) = IGolenetskii et a l. (2008 ^), , 

(l2008d ). (IT) ~ ^Golcn etskji et all (l2008dl} [12) - 

'Gplenctskii et al. (2008e), (IS) = lGolenetskh et all (I20q8g|'). (14 ) 
= IGolenetskii ct al .. (i2008hl'). (15 ) = IGolen etskii et al.l (120081 . 
(16) = iGolcnctskii etai.! (12008 ih. (17) = Fvon Kicnhn (2009^ 
(18) = Golcnctskii et al. (2009a'), (19) = Nakaaawa ct aD 
20091) ■ (20) = Golcnctskii ct al. (2009b.) ■ (21) = van dcr Hors^ 



20091 ) 
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Fig. 7. Correlation between PICsIT measured fluence and 
observed fluence from other instruments for 20 GRBs in 
our sample in three energy bands. Lines give the expected 
fluences derived from the conversion parameters in Table S) 



Table 4. Conversion coefficient from PICsIT counts to 
physical units in three energy bands. 



Energy K a 

(kcV) ( 10'^ erg cm~^ ct~^) 

260-364 09 03 

364-780 1.6 0.6 

780-2600 5.7 2.0 



bursts with higher fluences were detected. No significant 
dependence on the azimuthal angle (f> is evident. 

With the exception of GRB 061122, whose light curve 
is compromised by telemetry saturation, we compared the 
fluences observed by PICsIT to those derived in the same 
energy range for 20 GRBs with known spectral parameters. 
The ratio of the fluences in physical and instrumental units 
shows a strong dependence on the off-axis angle (see Fig.[6l 
upper panel), which can be partly compensated by a simple 
cos 8 factor, accounting for the geometrical projection ef- 
fect of the PICsIT detection plane (see Fig.|6l lower panel). 
The spread (in the worst case of a factor ~ 3) is caused 
by a more complicated energy- and direction-dependent 
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Fig. 8. Hard and soft color for 35 bursts in the PICsIT 
sample. HI = F2/F1 and H2 = F3/F2 are derived as 
the ratio of the physical fluences in 260-364 keV (Fl), 
364-780 keV (F2) and 780-2600 keV (F3). Diamonds refer 
to PICsIT instrumental quantities, and stars are obtained 
from documented spectral parameters. 

instrumental response as well as to the spectral parame- 
ters assumed for each single GRB. Neglecting the residual 
dispersion, the strong correlation shown in Fig. [7] between 
physical and instrumental fluences corrected for the geo- 
metrical projection can be assessed by a simple conversion 
coefficient: 



F{AE) = k{AE) X 



F, 



PICsIT 



(AE) 



COS0 



erg cm 



(1) 



computed as the weighted average of the ratio of the two 
quantities. The derived values of k in the three selected 
energy bands are given in Table HI with errors a estimated 
as the standard deviation. 

The rough calibration derived above allows us to com- 
pare the number of GRBs in our sample with the ex- 
pected burst rate based on results from previous exper- 
iments. From the peak flux cumulative distribution in 
the 260-364 keV energy band (Fig. we can estimate 
a completeness flux limit at ~600 ct/s, where the slope 
of the distribution flattens because of the decreasing effi- 
ciency to detect fainter GRBs at different off-axis angles. 
Assuming cos 0=0.5, and a Band spectrum with average 
valu es of the parameter s a= -1, 13— -2.3, and Epeak^SOO 
keV (jKaneko et al.ll2006|) , this count-rate limit corresponds 
to a flux of approximately 1 ph cm^^ s^^ in the 50-300 keV 
energy domain. From the BATSE LogN-LogP in this en ergy 
range (iKommers et al.ll2000l iFenimore fc Bloomlll995r ) the 
expected rate of GRBs with flux above the PICsIT limit 
is ~10 events/yr over a solid angle of An. This agrees rea- 
sonably well with the PICsIT detection rate of ~5.5 bursts 
per year (13 bursts in a net exposure time of ^850 days) 
in approximately half sky. 

Because most of the dependence on the burst incom- 
ing direction is accounted for by the detector area pro- 
jection coefficient, the ratio between fluences in different 
energy bands is almost direction-independent. We defined 
the quantities HI = i^2/Fl and H2 = F3/F2, where Fl, 
F2, and F3 are the fluences observed by PICsIT in the 
three selected energy bands, multiplied for the conversion 



Table 5. Burst variability: mean hardness ratio, x^, num- 
ber of bins and associated null-hypothesis probability. 



Burst 


Const 


x' 


N 




> xibs) 


H21 


060928 


1.6 


28.4 


14 


8- 


10"'' 


070829 


1.4 


15.0 


9 


6- 


10-2 


071003 


1.6 


22.9 


12 


2 • 


10-2 


080303 


1.0 


106.4 


23 







0803 19B 


1.0 


37.0 


21 


1 • 


10-2 


080615 


1.0 


17.5 


9 


2 • 


10-2 


080721 


1.3 


19.0 


13 


9- 


10-2 


080723B 


1.0 


38.0 


13 


2 • 


10-" 


H32 


060805B 


0.3 


13.8 


8 


5- 


10-^ 


060928 


0.4 


120.1 


13 







071003 


0.4 


18.5 


11 


5- 


10-2 


080303 


0.3 


46.4 


17 


8- 


10-^ 


080319B 


0.3 


44.6 


15 


5- 


10-^ 


080721 


0.4 


53.2 


10 


6- 


10-* 



coefficients in Tabled Figure [5] shows (red diamonds) the 
hardness ratios HI and H2 for the 35 bursts in our sample 
with non-null flucncc in the 780-2600 keV band, including 
also those events without localization. Star (blue) symbols 
mark the same quantities derived from spectral parameters 
of GRBs in Table H The region populated by PICsIT ob- 
servables is well confined within the expected values and 
characterizes the spectral colors of a typical burst observed 
by PICsIT. 

For those spectra described by a Band function we in- 
spected the dependence of instrumental hardness ratio on 
the spectral parameters, but we found no straightforward 
correlation of the soft color with the a parameter, nor for 
hard color with the f3 and Ep parameters. 

3.2. Spectral variability 

While a spectral characterization of a burst can be at- 
tempted only for events with known incoming direction, 
the hardness ratio evolution along the event is not depen- 
dent on the off-axis angle. From the background-subtracted 
count curves in the three selected energy bands 260-364 keV 
(soft), 364-780 keV (medium), and 780-2600 keV (hard), we 
computed the hardness ratios H21 (medium to soft) and 
H32 (hard to medium) along each burst. The hardness ra- 
tios were derived over the Tgg interval only for bins where 
counts are at least 2 • a above the background. The spectral 
variability was estimated by a simple test: 



X 



{HR, - {HR)f 



(2) 



where (HR) is the mean value over the N bins verifying 
the condition above, <7i is the error on the hardness ratio 
in the i— th time bin. When allowed by the signal-to-noise 
ratio, the hardness evolution for each burst is plotted in 
Figs. together with the event light curves. Note that 
a different binning time is usually adopted for light curve 
and HR plots. 

Table [5] reports the mean hardness ratio, the number 
of bins N, the value of and the associated null hypoth- 
esis probability of no-variability. Only bursts with a x^" 
probability higher than 90% and at least 8 bins are given. 
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Despite the strong Hmitation owing to the low SNR in 
short time bins, we identified nine events with significant 
variabihty {P > 90%), among which six are known GRBs 
and three are bursts without a documented position. The 
spectral variability in both hard to medium and medium to 
soft bands does not show a common evolution, as observed 
in different energy r egimes by other instruments (see e.g. 
iGuidorzi et aI1[201ll) . 

4. Conclusions 

We have investigated the capabilities of the IBIS/PICsIT 
instrument applied to the study of GRBs. We analyzed 
spectral timing light curves over about three years, pro- 
vided by PICsIT with energy resolution of eight channels 
and time sampling of ^^16 ms. Under stringent selection cri- 
teria, optimized for the search of long events, we identified 
39 bursts, most of which are documented GRBs or events 
detected by other telescopes. We showed that PICsIT suc- 
cessfully detects GRBs extending up to ~ 3 MeV at a rate 
of 5.5 events/year above the instrumental threshold of 
~ 600 ct/s in the 260-364 keV energy band. Bright events 
are observed even extremely off-axis, well beyond the nom- 
inal FoV. For all detected events we produced light curves 
with a time resolution down to 0.25 s. We used the spec- 
tral properties available in the literature of a subsample of 
GRBs to attempt a rough calibration between fluences in 
physical and instrumental units by separating the off-axis 
and energy dependence and computing an average calibra- 
tion coefficient in three selected energy bands. For suffi- 
ciently intense events the fine time and energy sampling of 
ST data allowed us to trace spectral evolution. We showed 
that it is possible to exploit the PICsIT archive for GRB 
spectral and temporal characterization in the hard X-ray / 
soft 7-ray regime. 
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Fig. 9. GRB 060805B background-subtracted light curve 
with 0.25 s binning time and hardness ratio with 0.5 s bin- 
ning time. Vertical (red) lines define the Tgo interval and 
horizontal (blue) lines show the average hardness ratios. 
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Fig. 10. Burst 060819 background-subtracted light curve 
with 0.25 s binning time and hardness ratio with 1 s binning 
time. See caption of Fig. IHlfor a description of the lines. 



Fig. 11. GRB 060901 background-subtracted light curve 
with 0.5 s binning time. 
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Fig. 12. Burst 060905 background-subtracted light curve 
with 10 s binning time and hardness ratio with 20 s binning 
time. See caption of Fig. [3] for a description of the lines. 
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Fig. 13. GRB 060928 - both episodes. Vertical (red) lines 
define the Tgo interval of the main episode. 
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Fig. 14. Second episode of GRB 060928: background- 
subtracted light curve with 0.25 s binning time and hard- 
ness ratio with 2 s binning time. See caption of Fig. [3] for a 
description of the lines. 



Fig. 15. Burst 061031 background-subtracted light curve 
with 1 s binning time and hardness ratio with 1 s binning 
time. See caption of Fig. for a description of the lines. 
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Fig. 16. Light curve of GRB 061122 in four energy chan- 
nels and with ^31 ms time resolution. Data gaps are rec- 
ognized in all energy bands. 
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Fig. 17. GRB 061222A background-subtracted light curve 
with 0.25 s binning time and hardness ratio with 0.5 s bin- 
ning time. See caption of Fig. [HI for a description of the 
hues. 
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Fig. 18. GRB 070207 background-subtracted light curve 
with 0.25 s binning time and hardness ratio with 0.5 s bin- 
ning time. See caption of Fig. |9] for a description of the 
lines. 



Fig. 19. Burst 070227B background-subtracted light curve 
with 1 s binning time. Vertical (red) lines define the Tgo 
interval. 
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Fig. 20. Burst 070326 background-subtracted light curve 
with 1 s binning time. Vertical (red) lines define the Tgo 
interval. 
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Fig. 21. Burst 070329 background-subtracted light curve 
with 1 s binning time and hardness ratio with 3 s binning 
time. See caption of Fig. [HI for a description of the hnes. 
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Fig. 23. Burst 070418 background-subtracted Ught curve 
with 1 s binning time and hardness ratio with 1 s binning 
time. See caption of Fig. for a description of the hnes. 




Fig. 22. Burst 070403 background-subtracted hght curve 
with 1 s binning time and hardness ratio with 2 s binning 
time. See caption of Fig. IHlfor a description of the hnes. 



Fig. 24. Burst 070429 background-subtracted hght curve 
with 1 s binning time and hardness ratio with 1 s binning 
time. See caption of Fig. [S] for a description of the hnes. 
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Fig. 25. Burst 070829 background-subtracted light curve 
with 1 s binning time and hardness ratio with 3 s binning 
time. See caption of Fig. [HI for a description of the hnes. 



Fig. 27. GRB 071003 background-subtracted light curve 
with 0.25 s binning time and hardness ratio with 1 s binning 
time. See caption of Fig. [S] for a description of the lines. 
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Fig. 26. Burst 070917 background-subtracted light curve 
with 0.25 s binning time and hardness ratio with 0.5 s bin- 
ning time. See caption of Fig. ^ for a description of the 
lines. 



Fig. 28. GRB 071006 background-subtracted light curve 
with 3 s binning time. Vertical (red) lines define the Tgg 
interval. 
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Fig. 29. Burst 071108 background-subtracted light curve 
with 0.25 s binning time and hardness ratio with 0.5 s bin- 
ning time. See caption of Fig. [HI for a description of the 
hnes. 
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Fig. 31. GRB 080204 background-subtracted hght curve 
with 0.25 s binning time and hardness ratio with 1 s binning 
time. See caption of Fig. for a description of the hnes. 
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Fig. 30. GRB 080122 background-subtracted Ught curve 
with 1 s binning time. Vertical (red) lines define the Tgo 
interval. 



Fig. 32. Burst 080303 background-subtracted light curve 
with 0.25 s binning time and hardness ratio with 1 s binning 
time. See caption of Fig. [S] for a description of the lines. 
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Fig. 33. GRB 080319B background-subtracted light curve 
with 1 s binning time and hardness ratio with 2 s binning 
time. See caption of Fig. [HI for a description of the hnes. 



Fig. 35. Burst 080408 background-subtracted hght curve 
with 1 s binning time. Vertical (red) lines define the Tgo 
interval. 
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Fig. 34. GRB 080328 background-subtracted light curve 
with 1 s binning time and hardness ratio with 2 s binning 
time. See caption of Fig. IHlfor a description of the lines. 



Fig. 36. GRB 0805 14B background-subtracted light curve 
with 0.25 s binning time and hardness ratio with 0.5 s bin- 
ning time. See caption of Fig. [3] for a description of the 
lines. 



V. Bianchin et al.: The first GRB survey of the IBIS/PICsIT archive, Online Material p 8 



080607 



08061 b 



15G 
100 
50 


200 
150 
100 
50 

-50 
100 

50 



Y:d 

0.8 
0.6 
0.4 

0.2 
0.0 





i 
¥ 


i \ 


1 






J 




\ J 






\ 












h 














ilf 








\ 










} ^ 



5 10 

Time (s) since trigger 




5 
Time (s) since trigger 



Fig. 37. GRB 080607 background-subtracted light curve 
with -.25 s binning time and hardness ratio with 0.5 s bin- 
ning time. See caption of Fig. [HI for a description of the 
hnes. 
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Fig. 39. Burst 080615 background-subtracted hght curve 
with 0.25 s binning time and hardness ratio with 1 s binning 
time. See caption of Fig. [S] for a description of the hnes. 
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Fig. 38. GRB 080613B background-subtracted Ught curve 
with 1 s binning time and hardness ratio with 1 s binning 
time. See caption of Fig. |9] for a description of the hnes. 



Fig. 40. GRB 080721 background-subtracted hght curve 
with 1 s binning time and hardness ratio with 1 s binning 
time. See caption of Fig. [S] for a description of the hnes. 
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Fig. 41. GRB 080723B background-subtracted light curve 
with 1 s binning time and hardness ratio with 1 s binning 
time. Data gaps are visible in all energy bands. See caption 
of Fig. [5] for a description of the lines. 



Fig. 43. Burst 080918 background-subtracted light curve 
with 0.25 s binning time and hardness ratio with 0.25 s 
binning time. See caption of Fig. [3] for a description of the 
lines. 
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Fig. 42. GRB 080817 background-subtracted light curve Fig. 44. GRB 090528B background-subtracted light curve 
with 1 s binning time and hardness ratio with 2 s binning with 2 s binning time. Vertical (red) lines define the Tgo 
time. See caption of Fig. |9]for a description of the lines. interval. 
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Fig. 45. GRB 090618A background-subtracted light curve 
with 1 s binning time. Vertical (red) lines define the Tgo 
interval. 



Fig. 47. GRB 090626A background-subtracted light curve 
with 1 s binning time and hardness ratio with 3 s binning 
time. See caption of Fig. for a description of the lines. 
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Fig. 46. GRB 090623 background-subtracted light curve 
with 2 s binning time and hardness ratio with 2 s binning 
time. See caption of Fig. IHlfor a description of the lines. 



Fig. 48. GRB 090809B background-subtracted light curve 
with 0.5 s binning time and hardness ratio with 0.5 s bin- 
ning time. See caption of Fig. [5] for a description of the 
lines. 



